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Abstract 

Direct numerical simulation (DNS), mostly used in fundamental turbu¬ 
lence research, is limited to low turbulent intensities due the current and fu¬ 
ture computer resources. Standard turbulence models, like RaNS (Reynolds 
averaged Navier-Stokes) and LES (Large Eddy Simulation), are applied to 
flows in engineering, but they miss small scale effects, which are frequently 
of importance, see e.g. the whole area of reactive flows, flows with apparent 
Prandtl or Schmidt number effects, or even wall bounded flows. A recent 
alternative to these standard approaches is the one-dimensional turbulence 
(ODT) model, which is limited to ID sub-domains. In two papers we will 
provide a generalized filter strategy, called XLES (extended LES), including 
a formal theory (part I) and one special approach in the XLES family of 
models, called ODTLES (in part II (see Glawe et ah (2015))). ODTLES 


uses an ODT sub-grid model to describe all turbulent scales not represented 
by XLES, which leaves the larger scales to be simulated in 3D. This allows 
a turbulence modeling approach with a 3D resolution mainly independent of 
the turbulent intensity. Thus ODTLES is able to compute highly turbulent 
flows in domains of moderate complexity affordably and including the full 
range of turbulent and diffusive scales. The convergence of XLES to DNS 
is shown and the unconventional XLES advection approach is investigated 
in basic numerical tests. In part II, highly turbulent channel and duct flow 
results are discussed and show the future potential of XLES and ODTLES. 
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1. Introduction 


Turbulence is a ubiquitous effect in the physical world, with major tech¬ 
nological and environmental impacts on human society and even human ex¬ 
istence. Turbulence inffuenced the density and associated gravity variations 
that led to the formation of present-day galaxies, stars, and planets. Without 
turbulent mixing, planetary atmospheric phenomena such as clouds, storms, 
and precipitation essential for life on Earth would be unimaginably different. 
Turbulence in the ocean has innumerable effects on oceanic biota, starting 
with the commingling of phytoplankton with needed nutrients. Due to its 
broad influence and baffling complexity, progress in fundamental and practi¬ 
cal understanding of turbulent mixing is exceptionally challenging, yet crucial 
for scientific advancement encompassing a wide class of problems in earth sci¬ 
ence, astrophysics, and engineering. For such problems, better understanding 
of turbulence interactions with buoyancy effects and chemical and thermody¬ 
namic processes is equally important. Due to its huge complexity, progress 
in understanding and prediction of turbulence is extremely challenging, yet 
crucial for scientific advancement in many disciplines. 

Analytic solutions of turbulent flows cannot be derived for real life appli¬ 
cations, e.g. with complex boundary conditions. Many meteorological and 
technical turbulent flows are investigated by numerical simulations. However 
the theoretical analysis of the governing equations is important to derive and 
improve simplifying models. 

There are several possible levels of turbulence modeling commonly ap¬ 
plied. The simplest model is no model; Direct Numerical Simulation (DNS) 
resolves all turbulent scales. Thus impacts of modeling and numerical errors 
are not influencing the investigated physics. DNS is widely used in funda¬ 
mental research, but limited to moderate Reynolds numbers Re (until now 
Re-j- 5200 for wall-bounded flows by Lee and Moser] (2014)) and Rayleigh 
numbers Ra due to the high computational effort, while many real-world 
flows have Rcr > 10® (following Smits and Marusic (2013)). 

In many industrial and some research applications Reynolds-Averaged 
Navier-Stokes equations (RaNS equations) are solved, because highly tur¬ 
bulent flows (e.g. Rcr > 10®) corresponding to realistic flows are computa¬ 
tionally feasible. RaNS describes the dynamics of time-averaged fields. The 
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influence of fluctuating terms is modeled. Various models are known (see e.g. 


Spalart and Allmaras (1992) or an overview by Pope (2000)). RaNS is used 


for steady states, but is generally not useful for computing time-accurate flow 
statistics. 

In recent years, Large Eddy Simulation (LES) has increasingly been used 
for industrial applied and fundamental turbulent flows. In LES spatially 
filtered equations are numerically solved while the unresolved sub-grid scale 
(SGS) terms are modeled e.g. by an eddy viscosity model (see e.g. 

eraL] ( pMj )). 


Germano 


Even in very fundamental highly turbulent flows LES need to resolve a 
wide range of scales including at least some portion of the inertial range of the 
turbulent cascade, which limits the achievable Reynolds numbers relative to 
the achievable Reynolds numbers in RaNS simulations. The parameterization 
of a certain range of small scales in RaNS and LES is especially problematic 
for multi-physics regimes such as buoyant and reacting flows because much 
of the complexity is thus relegated to the unresolved small scales. 


For further details on LES we refer to standard literature (e.g. by Sagaut 
( 2006| )). For an overview of turbulence properties and other model ap¬ 
proaches we refer to the work by Pope ( 2000| . 

Alternative model approaches e.g. the one-dimensional turbulence (ODT) 
model (see e.g. Kerstein (1999) and Kerstein et al. (2001)) describe the 


3D turbulence within a ID sub-domain including the full turbulent cascade, 
whereby the numerical representation of molecular diffusive effects becomes 


computationally feasible also in highly turbulent flows. Recently [Meiselbadi 
(2015) described wall bounded flows with Rcr < 6 x 10® using an adaptive 
ODT version (by Lignell et al. ( 2013[ )), which is clearly in the range of real- 
world applications, but limited to applications that are reasonable within a 
ID sub-domain. 

To beneht from the ability of ID models to describe highly turbulent 
flows, several approaches combine ID models with LES, e.g. LES-ODT (e.g. 
by Gao and Echekki (2008)), LES/LEM (by Menon and Kerstein ( 2011| )), 
and LEM3D (e.g. by Sannan et al. (2013)). Since ODT fully resolves e.g. 
molecular diffusion effects, it is fairly non-trivial to include it into an LES-like 
hlter approach. 

In this work we introduce an extended LES (XLES) hlter approach, which 
is tailored to include ID models as sub-grid models and can describe highly 
turbulent hows in a domain of moderate complexity, relevant e.g. in at¬ 
mospheric science. To achieve this, XLES solves 2D hltered equations on a 
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structured grid, maintaining one highly resolved Cartesian direction (this e.g. 
allows resolved molecular diffusion within an ODT sub-grid model). To de¬ 
rive a preferably general model, all Cartesian directions are treated equally: 
Three 2D hlters, each corresponding to one highly resolved Cartesian direc¬ 
tion, are applied independently to the governing equation. This leads to 
three coupled sets of 2D hltered equations, derived in section 

Thereby XLES only needs to resolve global 3D structures (e.g. the com¬ 
putational domain) independent of the turbulent cascade, while e.g. a one¬ 
dimensional turbulence model oriented in the respective highly resolved di¬ 
rections efficiently exploits the special symmetry of a 2D filter and represents 
the XLES microscale terms, corresponding to turbulent velocity scales not 
resolved by XLES. This ODTLES approach, investigated in detail in part II 
(Glawe et ah (2015)), is benehcial for flows with turbulent small scale effects 
playing an important role, e.g. in the combustion area, in buoyant stratihed 
flows and wall-bounded flows with high turbulent intensity e.g. with high 
Reynolds or Rayleigh numbers. 

This XLES based ODTLES model is an extended version of the ODTLES 


model, introduced and examined by Schmidt et ah (2008), Gonzalez-Juez 


et al. 

(2011 

), and 

Glawe et al. 

(2013 


In this work we distinguish between the expressions XLES to describe 
the XLES approach including an approximation or model for arising mi¬ 
croscale terms, ODTLES if these terms are in particular described by ODT 
and XLES-U (XLES unclosed), if microscale terms are neglected. 

In this work the XLES approach is introduced (in section]^ and XLES-U 
is verihed by numerical studies and applied to a turbulent channel flow (in 
section]^, followed by conclusions in section]^ In Part II (see |Glawe et al 


( 2015[ )) ODTLES is derived as one special approach in the XLES family 
of models and highly turbulent ODTLES channel and duct flow results are 
presented. 

Note that in this work no Einstein summation convention is used. 


2. Extended Large Eddy Simulation (XLES) 

We consider the turbulent flow to be described by the incompressible 
Navier-Stokes equations for a Newtonian fluid including the conservation of 
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( 1 ) 


mass in Eq. 0 and momentum in Eq. Q: 

3 

0 = ^ 


Hi “h ^ ^ dxjHj * Hi (2) 

i=i 

with the velocity Ui in Cartesian Xj-direction (z = {1,2, 3}), a constant kine¬ 
matic viscosity v, and the pressure gradient For simplicity we assume 
in this work the constant density to be po = 1 - 

Variables e.g. the velocities Ui and operators are assumed to be described 
within a continuum. Discrete variables are marked with the superscript d 
(e.g. discrete velocities: uf). 

Similar to LES, in XLES the velocity field is filtered. 3D filter functions 

commonly applied in LES are defined as tensor products of ID filter func¬ 

tions [Ik] (we are using an operator notation) in x^-direction {k = {1,2,3}; 
a filter operating on a continuum corresponds to a convolution). To derive 
macroscale (to be simulated) and microscale (to be modeled) terms the spa¬ 
tial scales are separated for modeling purpose using these ID filter operators: 

Ui = [ll 1] Ui (3) 

= [(/i -1- (1 — /i))(/2 -l- (1 — h)){h + (1 — ^ 3 ))] Ui 

= [{h + Si)(/2 + S2){l3 + -Ss)] Ui 

= [hhh + S1I2I3 + I1S2I3 + + -5152/3 -l- S1/2S3 -|- /1S2S3 -f- S1S2S3] Ui 

with the unity operator [1] and the ID small scale operator 

[sk] = [l-lk]- (4) 

The tensor product ansatz causes combinations of ID operators to be corn- 
mutable, e.g. S 1 I 2 I 3 = / 3 S 1 / 2 . 

In general the occurring spatial 3D scales (terms in the last row of Eq. 
(|^) are decomposed into those to be simulated (resolved) and those to be 
modeled (unresolved). We will refer to this as ‘filter separation’, because 
the commonly used term ‘scale separation’ implies the simulated scales to be 
large. In LES the resolved velocities correspond to the 3D large scale field 
[hhh]ui. Thus the LES filter separation corresponds to a 3D scale separation. 


0 = —d^^p -F 
Po 


dt-p^d: 

.7 = 1 


5 




Table 1: Comparison of filter separation ansatzes 


Model 

resolved scales 

unresolved scales (SGS) 

LES 

= [hklslui 

.yLES (s^/ 2/3 -)- I1S2I3 + 

ES1S2I3 + S1I2S3 + /1S2S3 + SlS2S3]Mi 

DNS 

[1 1 l]Ui 

0 

XLES 

Vj-i 

[hhh + S1I2I3 + I1S2I3 + hhsslui 

Ui = 

[S1S2/3 + S1/2S3 -|- I1S2S3 + SiS2S3]Mj 


In table [T] the filter separations are compared between LES, DNS and 
XLES. 

The LES resolved (large scale) velocities are and the 

unresolved (SGS) velocities = [1 — lil 2 h]ui- The acronym LES is used 
to differentiate between standard LES velocities and those dehned in XLES. 

In DNS all available scales are resolved numerically. 

In XLES the resolved scales are connected to 2D filtered fields. Without 
loss of generality (w.l.o.g.) we apply [/ 1 / 2 ] leading to: 


\l1l2\ui — [/1/2/3 + 


( 5 ) 


To treat all Cartesian directions equally, all possible 2D-filtered terms 
and [lil‘ 2 \ui are resolved numerically. 

This XLES macroscale (to be simulated) approach can be interpreted as 
the LES macroscale ([/i/ 2 ^ 3 ]) with additionally ID resolved small scale (RSS) 
terms in all Cartesian directions ([si/ 2 ^ 3 ] + [^ 152 ^ 3 ] + [^iG^s]). Both, the XLES 
resolved velocity scales (investigated in this work) and the unresolved scales 


(investigated in part II (Glawe et ah (2015))) impose special requirements 


on the numerical scheme and modeling; techniques known from LES cannot 
be applied one-to-one to XLES. 


2.1. XLES: Spatial Filtering 

In XLES, 2D filter operators -each corresponding to one highly resolved 
Gartesian direction- are applied to the governing equations, comparable to 
the 3D hltering in LES. To derive a hlter approach, independent of the cho¬ 
sen Gartesian direction, three 2D Liters, corresponding to three Gartesian 
directions, are applied. We will use a vector notation (called XLES vector 
notation), indicated by an underbar (e.g. uj. Each vector element contains 
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(a) LES grid: (b) Grid 1: (c) Grid 2: (d) Grid 3: 

~ [hl2h]ui- ^ [hlslui- ^ 2,1 ~ [hh]ui- uf i ~ [hl2]ui- 


Figure 1: In XLES the velocities are resolved using multiple XLES-grids 
illustrated in [Tbpd} 3D large scale properties, corresponding to a standard 


LES grid are for illustration represented with IVles = 4 cells per direction 
in Im In XLES the 3D large scale velocities are derived by ID hltering the 
XLES properties: = [lk]uk,i, k = {1,2,3}. The discrete XLES resolved 

small scale (‘RSS’) properties are represented for illustration by iV^ss = 16 


cells in lb Id 


one of the three 2D hltered velocity helds: 


[kh]ui = 0 

J/i/sjui/ \ 0 



( 6 ) 


with the 2D hlter matrix . 

The 2D hltered velocity helds (we refer to as the XLES vector) are 
discretized using three overlapping staggered XLES-grids with face-centered 
velocities and cell-centered pressure, illustrated in hgure [TbHTd] where each 
XLES-grid k = {1,2,3} discretizes one XLES vector element The size 
of the discrete ID hlter [If] corresponds to the large scale cell size (see 

hgure la), where the discrete ID hlter dehnition [If] corresponds to a box 
hlter 


_ 


[lf]Ui — 


^^LES 


Aj;LES 

k 


Ui dx'i.. 


(7) 


Here the discrete mesh corresponds to the discrete hlter size which is often 
called implicit hltering. 

The XLES hlter separation decomposes the full velocity held into three 
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parts (using the XLES vector notation): 


/ui\ /[l2h]ui\ /[hs2h + hhssluA /[‘SjnA 

Mi = Mi = [hl3]Ui + [sikh + hkS3]Ui + [5]Mi 

\Ui/ \[kl2]UiJ \[siy3 + hS2h]Ui/ \[‘5]tti/ 

= lAi + + [5]Mi : (8) 


1. ‘Directly Resolved’: 

The 2D hlter applied to the full velocity held leads to the XLES vec¬ 
tor: Mj = Ui- Each XLES-grid represents its own directly resolved 
velocity held distinct from the other XLES-grids. 

2. ‘Indirectly Resolved’: 

There are indirectly resolved small scale terms (‘directly resolved’ by 
another XLES-grid): Especially within non-linear advec- 

tion terms (in Eq. ([^), they determine the coupling between the XLES- 
grids. The coupling matrix C and the small scale matrix in Eq. 
^ are: 


C = 


/O 

1 

1 ^ 

\ i 

(Si 

0 

0 \ 

1 

0 

1 

and = 

0 

S2 

0 

ll 

1 

0 ) 

1 - 1 

^0 

0 



(9) 


The matrix dehnes the resolved small scale (RSS) velocities Mj 
(these are part of the XLES macroscale, contrary to LES): 

Ui = lki = iLi- = iki- ( 10 ) 

with the ID hlter matrix 

(h 0 0\ 

= 0 k 0 . ( 11 ) 

Vo 0 I 3 J 

In index notation: The term (w.l.o.g.) [/iS 2 / 3 ]Mi = [/i/sjMj —[/i/ 2 / 3 ]'Ui can 
be interpreted numerically, because M 2 ,i = [/i/ 3 ]ni is exclusively avail¬ 
able in XLES-grid 2 and the ID hltered XLES velocity held [l 2 ]u 2 i = 
[hkhl^i, corresponding to the LES velocity held (see hgureQ, is 
also available in XLES-grid 2. 


3. ‘Not Resolved’: 

The velocity scales [iSjwi = [S 1 S 2/3 + S 1 I 2 S 3 + / 1 S 2 S 3 + siS 2 S 3 ]Mi are not 
resolved in any XLES-grid (XLES sub-grid scale (SGS) or XLES mi¬ 
croscale). In part II (Glawe et al. ( 2015[ )) advection terms containing 
unresolved velocity scales are interpreted by the ODT model. 

We summarize all resolved velocity scales (directly and indirectly): 

% = [hkh + Slhh + hS2h + hhsslu, = (f ^ (12) 


A possible interpretation of the XLES macroscale is the numerical approx¬ 
imation of Ui instead of in LES (compare to table [^. [Ap¬ 


pendix 7 ^ shows that expressing in index notation automatically leads to 


three coupled XLES-grids (corresponding to Eq. (12)). 


2.2. Momentum Conservation 

Similar to LES, in XLES filtered equations are solved. Gontrary to LES 
the XLES equations are derived by applying a 2D hlter matrix to the 
momentum equation Eq. (j^ leading to (in XLES vector notation): 


0 =d^,P + 

i=i 


u. 

— I 




XLES 


(13) 


i=i 


1=1 


with the 2D filtered pressure 

t = = ([^ 2 / 3 ]^ [hh]p [hk]?)^ , 


(14) 


the XLES residual stress tensors and the Hadamard operator *, an 

entry-wise multiplication between XLES vectors and matrices. 

The three 2D filtered momentum equations, each describing three velocity 
components (hence the number of equations is tripled), coexist and are solved 
simultaneously. 

Gompared to LES advection terms the XLES advection terms 

(Wj^^ + hj) * + Ui) contain additional XLES-grid specihc 

resolved small scale (RSS) terms. 

The contribution associated with the XLES residual stress tensors 


XLES 

Lij 


{dxMj * Mj) - dx^u. * U. 


(15) 
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is captured through the use of some form of modeling or approximation. 


To investigate the residual stresses in more detail, a 2D decomposition (a 
modihed version of the 3D decomposition by Leonard (1975)) is performed, 
leading to: 


XLES 

—ij 


vXLES 


I ^XLES 
' —ij ' —ij 




(16) 


where: 


• XLES coupling tensor terms 


(uj »e »4, + »e 

( 17 ) 

involve ‘indirectly resolved’ terms. These stress terms couple the mo¬ 
mentum equations represented by different XLES-grids. A possible 


2 . 2.1 


In LES these 


approximation of is investigated in section 

terms are not simulated and therefore typically modeled or approxi¬ 
mated and contribute to 3D cross-stress and 3D SGS Reynolds stress 
terms. 


2D Leonard stresses 


£20 = {uj * Ui) - * uj + (f 


2 D vXLES 

■ij 




XLESn 


(18) 


describe the influence of the 2D hlter on the (‘directly’ and ‘indirectly’) 
resolved XLES advection terms (see section [2.2.2 for details). 


XLES cross-stress terms 


^xles 

—ij 


{uj *iki+ % * Mi) + {uj * + C_s 


(19) 


describe the interaction of resolved (‘directly’ and ‘indirectly’) and un¬ 
resolved XLES terms (see section 2.2.2|for details). 
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• 2D SGS Reynolds stresses 




( 20 ) 


describe the interaction of terms not resolved in XLES (see section 2.2.2 
for details). 


The XLES stress terms 


2 D /^XLES 


—ij 


and TZij 


XLES 


are investigated in more 
detail in section |2.2.2| and can be interpreted in the context of the ODT 
model, as shown in part II (see Glawe et al. (2015|)). 


A 2D decomposition following the idea of Germano (1986) is also possible 


within the XLES framework. Nevertheless the 2D decomposition following 


Leonard (1975) is snfficient to introdnce ODT into the XLES framework and 


therefore used here in preparation for part II (Glawe et al. (2015)). 


2.2.1. XLES: Resolved Advection Terms 

The coupling stress terms in Eq. (17) are decomposed into two 

terms (using from Eq. 


^xLEs ^ (a,»ft.: - ay® 4, (la^ »&)) 

( 21 ) 

with the matrix of ones 1 and where: 

1. The hrst term in Eq. (21) exclusively contains advection terms with 


advecting and advected velocities represented within the same XLES- 
grid. This can be interpreted as a linearization affecting the coupling, 
while the individual advection terms dx^Uk^Uk^i and in each 

XLES-grid k remain non-linear. We investigate w.l.o.g. one element of 


a transposed coupling vector [C_dxjUj * Uj) 
1 we hnd e.g. 


in detail: For XLES-grid 


dxj{[hh]uj [hh]ui) ( 22 ) 

with the velocities [lil^luj = U 2 ,j (j = {1,2,3}) only discretely repre¬ 
sented in XLES-grid 2. 

This transposed coupling vector can be rearranged: 

(£ dxjUj * Ui)'^ = * iki (23) 
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with the matrix 



( 24 ) 


Using Eq. (23) the same example in Eq. (22) leads to: 

(25) 

which is computed in XLES-grid 2 and then coupled to XLES-grid 1. 
Eq. (24) contains deconvolution operators which are realizable if 
only large scale information is present. This is the case in XLES: the 
velocities [lih]uj and \lih]ui in the example in Eq. (25) are large scale 
in the Xi-direction, which is the direction of the deconvolution. This 
implies that Eq. (23) is exact for a continuum. 

A discrete (numerical) approximation is provided by the algo¬ 

rithm shown in section 13.11 

Using the matrix 1} (including deconvolutions) the linearized coupling 
stress terms are: 




XLES 


= ij * Qp^. (uj * fl, - uy® * up®) 


( 26 ) 


2. The second term in Eq. (21) can be expanded as follows: 

dxj {lUj * liki - Uuj * hj) (27) 

+ UijUsp + U2,jUl,i + U2jU3,i + U^jUip + U3jU2p)- 


These terms contain interactions of small-scale velocities resolved in 
different XLES-grids. We neglect these non-linear coupling terms. For 
an intact energy cascade within the turbulent flow this assumption is 
reasonable because it implies that the velocities Ukj = Ukj — 
smaller than Ukj. 


The deconvolution within the coupling terms Eq. (26) is fundamentally 
different in XLES than in existing hltered treatments that do not have re¬ 
solved small scales. The reason is that the deconvolution in XLES is not 
intended to construct small-scale features that are otherwise non-existent, 
but rather, to modify a small-scale structure that already exists at the re¬ 
solved small scales, recognizing that this structure also has low-wavenumber 
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content. Indeed, the goal in principle is to modify appropriately that low- 
wavenumber content while preserving the high-wavenumber content to the 


greatest possible extent (see the example in Appendix C.3). 


By neglecting the non-linear coupling terms a XLES model error is intro¬ 
duced: 


“Axle's =4, (1% * ISi -1(% * m )) 


(28) 


In LES-U (unclosed LES), the LES-limit of XLES-U (see section 3.2[ ), 
the model error ct^les^ vanishe s an d thus comparing a convergence study of 


XLES-U and LES-U in section 


3.4) allows estimation of ct^les*- 


2.2.2. XLES-SGS: Leonard Stress, Cross-Stress and SGS Reynolds Stress 
Terms 

In this section the (to be modeled) stress tensors containing XLES mi¬ 
croscale (unresolved) terms -f- are 

investigated: 

• The Leonard stresses {uj *iLi + —dxj {ikj *iki + 

can in principle be calculated by explicit hltering of the directly and 
indirectly resolved XLES advection terms. 


The 2D Leonard stresses are expected to be small within a discrete 
formulation, since each 2D hlter size matches the corresponding XLES- 
grid size, similar to implicit hltering in LES. 


(and authors therein) report that explicit 
hltering the LES advection terms produces immoderate dissipation and 
neglecting the 3D Leonard stresses improves the outcome. 


Ferziger and Peric (1999 


According to these arguments we neglect 

Other attempts to model the 3D Leonard stresses are not transferred 
to the 2D Leonard stresses in this work. 


• The cross-stress and SGS Reynolds stress terms 


describe interactions including the XLES unresolved scales = [5]ui 
(see table [^. These terms can be modeled (e.g. by ODT). 

We decompose these stress terms: 


C 

—o 


XLES 


+ ^JLES ^j2D + |t „ £ ^ (29) 

+dxj (f ^ ) {Uj * *Ui + U^* Uj) 
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into: 

1. The terms dx^ + P * (uj * u^ + Uj * u^) 

depend on the resolved small scale velocities Uj exclusively avail¬ 
able in one XLES-grid. 

2. The terms dx^l^^ {uj * -I- + Uj) 

are independent of the XLES-grid (equal in all XLES-grids). If a 
modeling approach is simultaneously applied in different XLES- 
grids, a SGS-coupling is required to guarantee a consistent 3D 
large scale held. Ad hoc introduction of an additional coupling, 
e-g- + 1^ * C_^^) {uj * *Ui + Ui* Uj ), would lead 

to double counting of small scale terms available in all XLES-grids. 
Nevertheless the exact relation 


{Uj * - 1 - *Ui + Ui* Uj) 

= ^ + lj * Ql^^) (uj * *Ui + Ui* Uj) . 


(30) 


avoids double counting due to the factor | (proof of | + 1} * QJ^^) 

by insertion). 


We reformulate and summarize the XLES microscale terms in Eq. (29): 

= nr+ra r)M,, ( 31 ) 

with 

(«,■»“i + a,» ij + 1 (% »ly® + sy®» ii + ii»%)) , 

(32) 


Each vector element Aikjj contains the unresolved terms (SGS) in 
XLES-grid k which are coupled to the other XLES-grids by the SGS 
coupling terms ^ * C 
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Summarizing section 2.2.1 and 2.2.2 we can write the XLES momentum 
equations: 




^ ^ I spatial 

^.Uj * U-+ ^XLES 


i=i 

spatial 

discs 


i=i 


i=i 

^2 i= *=+ ij * At,j). 


(33) 


i=i 


The last line in Eq. (33) corresponds to the full coupling, including the SGS 
terms and the indirectly resolved terms. 

Additionally to the error term o^x^les (28) neglecting the Leonard 

stresses introduces the SGS error term: 


Spatial _ ^2D 

discs “ Aij • 


(34) 


Please note at this stage no concrete SGS model is introduced to approximate 




2.2.3. XLES: Mass Conservation 

In the incompressible flow regime, the hltered velocity helds need to be 
divergence free to ensure mass conservation. 

The 2D hlter matrix is applied to the mass conservation equation Eq. 
([^ (in the XLES vector notation): 

2=1 2=1 2=1 


which is decomposed into 3D large scale velocity fields and resolved 
small scale (RSS) velocities u^. 

A possible approach to solve Eq. (35) is to ensure mass conservation for 
both decomposed velocity fields and = s}^u{. 


1. The equation 0 = 

corresponds to the mass conservation in LES schemes. Standard dis¬ 
crete approaches applied in LES can be used. In our implementation a 
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pressure Poisson equation is solved. This leads to a large scale pressure 


= whose gradient enforces a divergence free velocity 


field 

field by solving = 0 (for details we refer to stan¬ 

dard textbooks, e.g. [Ferziger and Peric (1999)). 

The equation 0 = 

corresponding to the RSS velocity fields, is discretely fulfilled without 
additional effort under some conditions. 

In [Appendix B the XLES mass conservation is derived in detail and 
three conditions are identified to ensure the RSS velocities to be diver¬ 
gence free: 


Consistency condition: is equal in all XLES-grids k which is 

valid due to coupling (see section 3.2 for details). 


• A divergence free 3D large scale velocity field which is valid after 
the standard pressure projection. 


• A discrete ID box filter defined in Eq. ([^ is used. Here we intro¬ 
duce a coarse-grained (large scale) and staggered control volume 
of the size (this cell size equals the filter size). Please note 

this filter definition is only used explicitly to create = [lk]uk,i 
within the coupling and the mass conservation. The 2D filter 
is not applied explicitly to any term. 


In summary the discrete mass conservation is assured by a standard 3D 
approach (large scale) with CI(Xles) cells, if a box filter and a staggered grid 
is used. 

In XLES the velocity fields, including consistent 3D large scale infor¬ 
mations, are discretely interpreted in three XLES-grids simultaneously. For 
such a system an existing redundancy of the discrete velocity fields can be 
exploited by deriving the velocity components ^ within one 3D large scale 
cell due to a direct solution, w.l.o.g. in XLES-grid 1: 


u 


1,1 




+ Xi = u 


_ —d, LES 


1 


Ax^es 


(36) 




A:cLES 


+X1 


5^,2 uf 2 dai'i 




A„LES 


+x\ 


gdx) 


for Xi < Ax^^®. Eq. (36) is a semi discrete interpretation of Gauss’s theorem 
for a divergence free velocity field tailored for velocities in XLES-grid 1 (see 
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figure lb) within one 3D large scale cell (see figure la) on a staggered face 
centered grid (no interpolation necessary). 


Since w.l.o.g. is specified due to Eq. (36), the momentum equa¬ 


tion Eq. (33) in XLES-grid 1 only needs to be solved for i = {2,3}. In 


consequence 6 momentum equations are dynamically solved (2 velocity com¬ 
ponents in each of 3 XLES-grids), while 3 velocity components can be derived 


by Eq. (36). 


These reduced momentum equations can be expressed by multiplying a 
Kronecker delta operator matrix 


'1-5 


0 


0 


t-5,= 


1* ^ \ I n if h = ? 

0 1 -^ 2 . 0 I ,with(l- 4 .) = <l; , ( 37 ) 

0 0 1 - 53V 


1 , else 


giving: 


0 =(| - 4 ) ( + ldt-u^dl.\u. + ^d. 


■ u. * u. 

j—] —t 


(38) 


i=i 




+(!-§ 


I + (1 - ^)a" 


0=1 


+4 - § E (f + 4 - fi) E k * . 

t=i i=i 

Please note the Kronecker delta operator matrix is also applied to the unre¬ 
solved terms and therefore can be exploited within a model approach. 
Additionally the factor ^ within Eq. (32) can be replaced by 


Since in XLES no 3D small scale velocity held is dehned, small scale 
pressure effects can only be captured by the modeling approach (additionally 
to the modeling of the SGS advection terms Alp)- 


3. Properties of the XLES Approach 

3.1. XLES: Deconvolution 


In sections 2.2.1 and 2.2.2 coupling terms between the XLES-grids are 
introduced. A discrete approximation of the transposed coupling matrix 
(G/^°)^ = is required to numerically represent these coupling terms. 
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has to be interpreted. 


Within the matrix (see Eq. (24)) an discrete deconvolution operator 


[I 

In general a velocity held u (indices are skipped in this section and all 
variables are assumed to be discrete representations) cannot be reconstructed 
exactly 7 ^ 1 ), because by hltering (e.g. = [l]u) information get 

lost. This is not recoverable, unless only large scale information is present 
in the full spectrum, which is fortunately the case in XLES, as shown by an 


example in Eq. (25). 


The deconvolution within XLES is not intended to construct small-scale 
features, but rather, to modify the low-wavenumber content of an existing 
fully resolved property. This modihcation must meet several requirements: 

• Integral constraints are imposed to satisfy conservation laws at the level 
of the individual 3D coarse-grained cells. This means the box hltered 
large scale held needs to be preserved by the deconvolution, e.g. to 
enable the consistency preservation of the 3D velocity held required for 
the mass conservation and the LES limit of the XLES model (in section 


3.2): 


['ir‘i = 1- 


(39) 


The deconvolution must in some sense modify the large scale structure 
along the resolved direction while in some sense preserving the small- 
scale structure. In this sense this deconvolution is a reconstruction of 
existing structure rather than a construction of something that does 
not otherwise exist. For this reason the term reconstruction has been 


used to describe it (Schmidt et ah (2008) and McDermott (2005)) 


In principle the deconvolution approach is not restricted to one particular 
hlter dehnition (e.g. spectral hlter, box hlter, Gaussian hlter), but the mass 
conservation is greatly simplihed if a box hlter is used, as shown in [Appendix 


In this section we introduce a discrete approximation of the deconvolution 
operator [/“^] for the purpose of computing a highly resolved velocity held 
Ua ~ (a indicates an approximation by discretely interpreting [l~^]) 

from a box hltered held = [l]u. 


Schmidt et ah (2008) introduces a recursive and very fast algorithm ap- 
-ii^yLES order accuracy without changing the 


proximating Ua ~ [/ 
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Figure 2: A representative instantaneous streamwise velocity profile u (solid) 
for a turbulent channel {Rcr = 1020) with A^les = 16 and A^rss = 1024 is 
shown. The large scale velocity = \l]u is represented with points. The 
deconvolution using the original algorithm by Schmidt et ah 


(2008) is dashed-dotted and the limited algorithm is dashed. For both algo¬ 
rithms (limited and unlimited) the reconstructed held Ua preserves the box- 
hltered values [/]«„ = but without the limiter unphysical overshoots 

occur near the walls. 


box hltered held {\l]ua = = [l]u , see Eq 

tails, we refer to [Schmidt et al. (2008), 


(39)). For algorithmic de- 


Gonzalez-Juez et al. (2011), and 


McDermott (2005). 


The algorithm produces artihcial local extrema in areas of monotone large 
scale helds with large gradients, e.g. occurring near a wall as reported by 


McDermott (2005). Similar ehects are known from deconvolution approaches 


within the helds of image and signal processing. We expand this algorithm by 
applying the linear slope limiter by Burbeau et al. (2001) after each recursive 
step in the area of monotone velocities. Hereby monotonicity is dehned 
locally for one cell by including its neighbors. 

The properties of the original algorithm and the limited one are demon- 
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strated for a realistic instantaneous streamwise velocity profile u for a turbu¬ 
lent channel case with high gradients at the walls, shown in Figure]^ (please 
note in XLES only large scale information is deconvolved). 

The numerical deconvolution error due to the algorithm described can be 
written formally: 




patial ={IUQ- * Q [u- *U.- u] 


*Ui 


+ f ""Mp) • (40) 


There is no spectrally sharp way to implement this type of deconvolution, 
Various non-equivalent procedures are possible in principle. [Schmidt et al. 


(2008) describe a particular approach that is adopted here with a technical 


modification that improves its behavior near walls. To summarize, although 
the inverse operator arises as a natural and necessary consequence of 
the XLES ansatz, the ansatz per se does not uniquely define its meaning nor 
guarantee that it can be specified in a way that is free of unintended artifacts. 


3.2. XLES: Consistency Preservation and ‘LES limit’ 


Both, the XLES mass conservation in section [2.2.3| and the coupling terms 
in section 2.2.1 and 2.2. 2| assume the 3D large scale velocity field to be consis¬ 
tent, meaning each XLES-grid k contains the identical 3D large scale velocity 


field: 




= u 


;rrLES _ 






In the XLES-vector notation this condition can be written as: 





(41) 


with the ID filter matrix defined in Eq. (0. This means consistency pre¬ 
serving terms should be independent of the XLES-grid k. 

We assume that the initial conditions are consistent (this can easily be 
achieved) and need to prove a consistency preserving XLES-advancement 
including advective and diffusive terms and the sub-grid modeling. 

The ‘directly resolved’ advection terms themselves {dx^Uj *Uj) are violat¬ 
ing the consistency condition, but by including the corresponding coupling 
(indirectly resolved terms): dx^Uj *Uj + the consistency is preserved. 

The same is valid for microscale terms and their couplings: * 
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Consistent 3D velocity field 


O 



Mass conservation! 

_I 



Coupling 


Figure 3: On the one hand the consistency preservation is invoke before the 
advection step to define the coupling terms. On the other hand the coupling 
guarantees the consistency to be still preserved after the advection step. 


By rearranging the XLES advection terms and their corresponding cou¬ 
plings, we find: 


Uj *Ui+ Kij 


XLES _ - , — 


= Uj * Ui 


(42) 


with = [lil 2 h + Sil 2 h + hs 2 h + which is equal in all XLES-grids 

k. 

The same holds for the XLES sub-grid terms and their corre¬ 
sponding couplings (proof by insertion). 

In consequence the coupling terms are essential for a consistent 3D veloc¬ 
ity field within XLES, as illustrated in figure 

In the implementation of XLES we check that the 3D large scale velocity 
field is still consistent every 100 timesteps to avoid numerical errors violat¬ 
ing consistency preservation. All computations show the consistency to be 
limited by the floating point accuracy. 

In the limit A^rss —t A^les the 2D hltered XLES-U equations collapse 
to the 3D hltered LES-U equations in each XLES-grid, because the RSS 
velocities vanish: iii = Ui — = 0 

For Ales Arss(> Adns) in each direction, all velocity scales are re¬ 
solved by the 3D grid and XLES-U (and LES-U) converges to DNS. 

Within the ODT model, applied as an SGM in part II (Glawe et ah 


( 2015| )), all turbulent events are suppressed in the DNS-limit. Thus DNS is 
also a distinguished limit of ODTLES. 


21 


















3.3. XLES: Coupled Advection Scheme 

The characteristic shape of the XLES-grids (hgure lb-Id) is considered 
when choosing the nnmerical schemes to be implemented, e.g. an implicit 
time discretization in the highly resolved direction is applied, while explicit 
time schemes are nsed to advect large scale properties. Therefore different 
nnmerical advection schemes are mixed dne to the XLES conpling terms, be- 
canse the same property, represented by several XLES-grids is simnltaneously 
advanced by different nnmerical schemes. Additionally the conpling reqnires 
a deconvolntion fnnction, which interacts with these nnmerical schemes. 

In this section we present the nnmerical properties of this conpled advec¬ 
tion scheme. For simplicity the 3D large scale velocity helds are resolved 
by Ales cells in each direction. The resolved small scale (RSS) properties 
are discretized with Arss cells in all XLES-grids. 

The conpled advection scheme nnmerically approximates the eqnation 


0 =(1 - dtu 


j=i 


■ U . *u 

j—j 


i=i 


[Uj 


*Ui 




(43) 


(according to Eq. (38); SGS terms are not considered). 

Fnndamental properties of the conpled advection scheme, e.g. the nn¬ 
merical dissipation and dispersion, can be demonstrated by solving a one 
dimensional linear advection problem. Then a constant wave speed cj re¬ 
places the advecting velocities in Eq. (43): ii^ = = Cj (ill). In 

this section the linearized conpled advection scheme is investigated, while 
the fnll non-linear coupled advection scheme is studied in section 

W.l.o.g. the advected velocity is represented by the two staggered 
XLES-grids containing U 2 ,i and its^i respectively while ui^i is evaluated using 


Eq. (36). 


Two situations occur (see hgure |^, as discussed here by a concrete ex¬ 
ample (generalization to other indices is trivial): 

• Case 1\ i = i l in Eq. (43): Advection into Xi-direction; In both 


XLES-grids 2 and 3 the advection is coarsely resolved with Ales grid 

cells. 


Case 2: 


2 j ^ 


. W.l.o.g. ^ 
I = \ 


in Eq. (43): Advection into direction 


X 2 ] In XLES-grid 2 the advection is resolved with Arss cells, respective 
with Ales cells in XLES-grid 3. 
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(a) Case 1: j = z in Eq. (43). Coarsely (b) Case 2: j / z in Eq. (43). Highly 
resolved advection in XLES-grid 3 (left) (coarsely) resolved advection in XLES- 
and XLES-grid 2 (right). grid 3 (left) (XLES-grid 2 (right)). 


Figure 4: XLES requires a coupled advection scheme. The advected velocity 
hj (bold arrows) is advected with the constant wave speed Cj (dashed arrows) 
and is resolved in two XLES-grids. For the two relevant cases the discrete 
staggered grids within one 3D large scale cell (compare to £gure[^ are shown 
in (a) and (b). 


To spatially discretize the advection terms, a central difference method 
(CDM) on a staggered grid is used. 

Two different time discretizations are deployed: 


The highly resolved properties (using N^ss cells) are discretized using 


a standard 2nd order Crank-Nicolson (CN) scheme in time (see Crank 


and Nicolson (1996)). We will refer to this fully discrete scheme as 


CN-CDM. 


The coarse grained resolved properties (using IVles cells) are discretized 


using a 3-stage 3rd order TVD Runge-Kutta (RK3) scheme (see Spiteri 


and Ruuth (2002)). We will refer to this fully discrete scheme as RK3- 


CDM 


The numerical properties of the numerical schemes without coupling are: 

• The advection scheme CN-CDM is stable, dissipation free and has low 
dispersion. 

• The RK3-CDM scheme is found to be stable, producing little dispersion 
and dissipation. 

These properties transfer to the coupled advection schemes (in Cases 1 and 

2 ). 

Based on the expressions RK3-CDM and CN-CDM the coupled advection 
scheme in Case 1 (Case 2) is called RK3-RK3-CDM (CN-RK3-CDM). 
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To illustrate the coupling procedure w.l.o.g. for the CN-RK3-CDM scheme, 


we linearize and rewrite Eq. (43) in a semi-discrete form (in index notation). 
Hereby is discretized in XLES-grid 3 in Eq. (44) and ■U 2 ,i in XLES-grid 


2 in Eq. (45): 


{At)u3^l = / C2d^^U3^i dt 

JrK3 


(44) 


+ [^3 ] [^2J \ I '^■2^X2 

'CN 

{At)u2,l = / C2dx2U2,l dt 

JCN 


029x2^2,1 dt 


C2dx2U^^ dt 


lRK3 


(45) 


+ [^2 ][y( / 028:^2^3,1 dt- I C 29 „ 2 M 3 i“dt 

'RK3 JrK3 


Here the numerical time discretizations of the advection terms are indicated 
(CN and RK3). The discrete coupling operators [Ik] and (see Eq. ([^ and 
the algorithm described in section 3.1) are discretized by an explicit Euler 


scheme in time. Also the non-linear advection part which is not considered 
in the linear advection problem is discretized by an explicit Euler scheme in 
time. 

The properties of the fully discrete linear coupled advection schemes are: 

• The coupled advection scheme converges to the analytical solution for 
^LES —t cxD as investigated in [Appendix C.2| where the theoretical pre¬ 


diction e.g. by Tsai et ah (2002) is reproduced. This property is re¬ 


quired to numerically realize the LES limit and DNS limit of XLES 
(section |3.2[ ). 

The coupled advection scheme including the deconvolution (see section 


3.1) reproduces well dehned large scale behavior, even by simultane¬ 


ously transporting small scale properties, as investigated in [Appendix 

Q 

The highly resolved CN-CDM scheme (using A^rss cells) contributes to 
the coupled CN-RK3-CDM scheme (Case 2) and additionally increases the 
numerical accuracy of the coupled scheme. This is neither required for a well 
dehned and converging XLES scheme nor the LES limit (or DNS limit) of 
XLES. As investigated in [Appendix C.2 the highly resolved advection terms 
increase the coupled numerical accuracy up to a resolution ratio A"rss/A"les ^ 
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10. With IVrss ^ IOIVles the overall numerical error is dominated by the 
coarsely resolved RK3-CDM scheme. 

The coupled numerical schemes, simultaneously discretized by multiple 
XLES-grids, are found to be appropriate for linear advection. Important 
numerical properties, like stability and dissipative behavior, transfer from 
the well known underlying numerical discretizations to the coupled advection 
scheme. 

There are other time schemes, that possibly can be adapted to the special 
requirements of the XLES time advancement, e.g.: 


An IMEX (implicit/explicit) time scheme (see e.g. Cavaglieri and Bew- 


ley (2015) and references cited therein) applied to XLES can in principle 
lead to a high-order time integration for all terms (including non-linear 
advection and coupling), but to include the coupling terms into such 
a high order time scheme requires identical coefficients for the implicit 
and explicit integration terms. 


Adapting split-explicit schemes (see e.g. Gadd (1978)) to the XLES 
coupled advection can especially decrease the dispersive effects that 


arise. 


Large time step wave propagation schemes based on the work by LeV- 


eque (1985) can perhaps improve the numerical properties of the ap¬ 


plied CN scheme within the XLES advection scheme. 


For the sake of completeness the diffusion terms are discretized by a 
central difference scheme in space with a hrst order explicit (implicit) Euler 
scheme in the coarse (hne) resolved XLES-grid direction. 


3.4- XLES: Convergence of the XLES Approach 

To verify the XLES model, we conduct a convergence study of XLES-U 
(unclosed XLES) for a fully developed turbulent flow, including the full diver¬ 
sity of non-linear advective effects: A turbulent channel flow with a friction 


Reynolds number Rcr = 395 is computed. DNS results by Kawamura et ah 


(1999) (online available: Kawamura (2014)) are compared to the XLES-U 
outcome. 

The DNS is resolved with N^ns = 192 non-equidistant cells in the hor¬ 
izontal direction (between the walls). For the spatial discretization a sec¬ 
ond order central difference scheme is used. The time is discretized using 
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a Crank-Nicolson scheme for the wall-normal non-linear terms and a second 
order Adams-Bashforth scheme for other terms (results for higher resolutions 
and higher order schemes are also available online). This numerical scheme 
is comparable to the XLES-U numerical scheme (see section 3.3). 


For the XLES-U convergence study, the number of equidistant 3D large 
scale cells Ales is increased using the values Ales = {16,32,64} while 
Arss = 512 is kept constant. 

The time step size is limited by: 


At = CFL min —A: = {1, 2, 3} 

\ Uk,i J 


(46) 


with a Courant-Friedrichs-Lewy number CFL = 0.45 and the small scale cell 


size 


To produce statistically signihcant results, the flow is averaged for tave > 
25 non-dimensional time units (compared to tave = 20 for DNS) after reaching 
a statistically steady state. 

Additionally equidistant discretized LES-U (LES limit of XLES-U with 


Arss = Ales, see section 3.2) channel flow results are compared to XLES-U. 


The latter resolves selected advective and diffusive small scale effects even 
without using a SGM. Those can be identihed by comparing the XLES-U 
and LES-U channel flow results. 

Figure illustrates the results: The mean velocity prohles computed by 
XLES-U (see hgure 5b) and LES-U (see hgure [5a| are compared to DNS. 

Additionally the streamwise and spanwise velocity RMSs (see hgure 
and the budget terms of the turbulent kinetic energy (see hgure [Sd] - [5f| ) are 
shown for XLES-U, LES-U, and DNS. 

Both XLES-U and LES-U show convergence towards the DNS results with 


increasing 3D resolution. Following arguments in section [372] this implies the 
XLES error term ^^xles small. Additionally the XLES non-linear 

terms (and its LES limit) are sufficiently represented and converging, even 


if the under-resolved velocity prohle (hgure 5a and 5b) indicates noticeable 
numerical dissipation. 

XLES-U is able to accurately represent dihusive ehects, e.g. the laminar 
sublayer near the walls, independent of the 3D large scale resolution and 
reproduces basic advective ehects including turbulence with Ales ^ 32 3D 
cells for Rcr = 395. Additionally the turbulence statistics proht from the 
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shifted with increasing shifted with increasing 



(c) Streamwise (ttRMs) and spanwise 
(^^RMs) velocity RMS. 


LES i 

XLES 

« 


"‘■v 


s 



-100 0 100 


y-^ 

(e) Viscous transport of the turbulent ki¬ 
netic energy (tv). 



(d) Production (Prod) and Dissipation 
(Diss) of the turbulent kinetic energy. 



(f) Advective transport of the turbulent 
kinetic energy (ta). 


Figure 5: Turbulent channel flow results for DNS (small crosses), LES-U, and 
XLES-U with A"les = 16 (dash-dotted), iVLEs = 32 (dashed), and Ales = 
64 (solid). The XLES-U small scales are resolved using Arss = 512 cells. 
XLES-U mean and statistical flow properties are based on the velocity field 
perpendicular to the channel walls U 2 i (in the LES limit: U 2 i = 
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additional small scale effects represented by XLES compared to LES-U even 
with low 3D resolutions. 


In part II (Glawe et al. (2015)) ODT is demonstrated to represent the tur¬ 


bulent effects not resolved by XLES-U, even with very coarse 3D resolutions 
(e.g. A^les = 16). 


4. Conclusions 

In this work XLES, an extended LES model, is introduced. This approach 
is intended as a basis for a new class of turbulence models, tailored to describe 
a (XLES specific) macrostructure of highly turbulent flows including fully 
resolved molecular diffusion in domains of moderate complexity. 

An innovative XLES filter strategy is derived, including a numerical rep¬ 
resentation of the resulting 2D filtered XLES equations by three coupled 
XLES-grids. This coupling requires a high order deconvolution function for 
ID box filtered fields, wherefore a stabilized version of the algorithm intro¬ 


duced by Schmidt et al. (2008) is applied 


The 3D hltered LES-U (unclosed LES) equations and DNS are distin¬ 
guished limits of XLES-U (unclosed XLES). 

We verify the XLES-U equations: 

• By analyzing the coupled (linear) advection scheme occurring in XLES. 

• By performing a convergence study: The turbulent channel flow with 
Rcr = 395 is compared to DNS and LES-U. Compared to LES-U low 
3D resolutions are required in XLES-U to reproduce fundamental flow 
statistics, especially in the near wall region. E.g. the laminar sublayer 
is fully represented independently of the 3D resolution. 

An XLES closure can be provided by the ODT model, which is able to 
capture the full turbulent cascade of highly turbulent flows, as introduced in 
part II (Glawe et al. ( 2015[ )). Recently Meiselbach (2015) introduces chan¬ 
nel flow results up to Rcr = 600000 using adaptive ODT by [Lignell et al. 


(2013). This ODTLES model enables the representation of highly turbulent 


flows in domains of moderate complexity with coarse 3D resolutions basically 
independent of turbulent scales and at once represents diffusive and (turbu¬ 
lent) advective effects down to the molecular level using ODT. Hereby the 
3D resolution needs to be sufficient to resolve the domain and possible sec¬ 
ondary effects. E.g. due to secondary instabilities within a duct flow, the 3D 




















resolution indirectly depends on the Reynolds number (investigated in part 
II by Glawe et al. (2015)). The XLES model is especially advantageous for 
turbulent flows with crucial small scale effects, e.g. with important Prandtl 
or Schmidt number effects, buoyant stratihcation effects, or in the held of 
combustion. 

In XLES a potential bottleneck in highly resolved incompressible simula¬ 
tions, the pressure handling, is only computed on a very coarse 3D grid and 
small scale pressure effects might be modeled. Closed XLES (e.g. ODTLES) 
is a promising and highly parallelizable approach to investigate fundamental 
atmospheric turbulent hows. 
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Appendix A. XLES: Vector Notation 

We introduce an alternative to the vector notation in section |^by writing 
the XLES velocity scales in index notation. The XLES vector includes the 
XLES vector elements Uk,i represented in XLES-grids k, k = {1, 2, 3}. 

In index notation the resolved small scale (RSS) terms Si/ 2 / 3 , / 1 S 2 / 3 , and 
/ 1 / 2 S 3 can be expressed by one term: Sklk®ih® 2 , where {k, fc © 1 , fc © 2 } is a 
positive permutation of {1,2,3} (the ID hlter operators are commutable). 
The operator © denotes the positive permutation: 

q (B X = {{q + X — 1) mod 3) + 2 (A.l) 

for q = (1, 2, 3} and x = (1, 2} (valid in three dimensions). 

The XLES resolved velocity scales (Eq. (|I^) are written in index notation 
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(using the operator ©): 


Ui — [/l/2^3 + S 1 I 2 I 3 + hS2h + ^1^2'S3] Ui 


kkh + E 


902 


9=1 


Ui 


(A.2) 


E(v g©l^ij©2 + Sg/g©l/g©2) — 2/1/2/3 
.9=1 


Ui . 


The 3D large scale operator [/i/ 2 ^ 3 ] in the last row of Eq. (A.2) can be 
expressed in terms of [/q/g©i/g ©2 + Sg/q©i/g© 2 ]: 


‘^kkk ^ ^ (1 ^qk){J‘qk®^k®‘^') 

9=1 

3 

'y ^ (1 Sqkjiikk k ^q^g)^g©l^g©2) 

9=1 

3 

y ^ (1 Sqk)kikk®^k®^ A (1 Iq) lq®llq®2) 


(A.3) 


9=1 


with the Kronecker delta operator dehned in Eq. (37). 

Here an additional arbitrary index k = {1,2,3} is introduced. Eq. (A.3) 
is satished for k = 1, k = 2, and k = 3. This index k spans the XLES vector. 

Insertion of the 3D large scale operator (Eq. (A.3)) into the XLES re¬ 
solved velocity scales (Eq. (A.2)) leads to: 

3 




y ^ (t kq ^g^gfc)) kqk®^k®‘^ Sqk®^k®‘^ 

L9=1 


Ui, k = {1,2,3}. 

(A.4) 


This XLES resolved velocity Uk,i reproduces exactly the XLES resolved ve¬ 
locity scales in Eq. (12) including the XLES coupling XLES is 

interpreted as approach hltering the Navier-Stokes equations by applying 
the operator [/i/ 2^3 + sikh -l- hs 2 h + ^ 1 ^ 253 ] compared to [/i/ 2 ^ 3 ] in LES (see 
table [^. 
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Appendix B. XLES Mass Conservation: Resolved Small Scales 


Divergence-free XLES velocity fields are guaranteed, if the 3D large 
scale velocity field = [lk]uk^i fulfills three conditions (see section 2.2.3): 

• Con.l: is consistent 

• Con.2: is divergence free, enforced by 3D standard approach. 

• Con.3: [Ik] is a discrete ID box hlter (dehned in Eq. ([^) in x^-direction. 

The XLES velocity helds u^ are divergence free as proved within one 3D 
large scale cell of the size Axk in x^-direction (operator © is dehned in Eq. 

(M): 


0 = 


3 

j=l 


+ Uk,j) dxk = 


- 3 

'> —^ j=i 


i'^kjj dxk 


(B.l) 


dxk © 


^^kS)l^k,k(Bl dXk I ^a:j.02^fc,A:©2 dXk 


^k,k 


Axi 


Cra.3 lES 


2 
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for each XLES-grid k = {1,2,3}. Here a staggered grid is used, leading to 
Uk,k (^f^) without additional interpolation. 

This in particular means the resolved small scale (RSS) velocity helds 
are divergence free by construction. 


Appendix C. XLES: Numerical Advection Scheme 

The XLES equations are discretized using multiple, coupled XLES-grids. 
This discretization includes a coupled advection scheme, possessing uncertain 
numerical properties. 
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The coupled advection scheme is investigated and verihed by computing 
a linear wave propagation. The initial condition 

Ml = sin(47ra;i) (C.l) 


and periodic boundary conditions are used. The number of cells for the 
resolved small scale helds A^rss and the 3D large scale held A^les are varied. 
The wave is propagated for 10 wavelengths {t = 5 L/cj; L = 1; constant wave 
speed Cj). 

All calculations are done with a Courant-Friedrichs-Lewy number CFL = 
0.45 limiting the time step size 


At = CFL min 

i,k 




,i,k = {1,2,3}. 


(C.2) 


Please note the non-linear problem in section 3.4 is limited to time steps 


based on instead of which is investigated in part II (Glawe 


et al. (2015)) in detail. 


Appendix C.l. 


Coupled Advection Schemes: RK3-RK3-CDM and CN-RK3- 
CDM 


Zha and Lingamgunta (2003) showed that RK2-CDM (central difference 


in space and 2nd order Runge Kutta in time) is unstable, while RK4-CDM 
(4th order Runge Kutta in time) is stable and fully dissipation free. 

We found the coupled advection schemes RK3-RK3-CDM (Case 1 in sec¬ 
tion 3.3) and CN-RK3-CDM (Case 2 in section 3.3) to be stable and slightly 


dispersive (see hgure [06 ). Additionally RK3-RK3-CDM is little dissipative. 
The difference between RK3-RK3-CDM and CN-RK3-CDM is not signih- 
cant. With increasing 3D resolution (Ales) fhe dispersive behavior decreases 
(convergence). 


Appendix C.2. Coupled Advection Scheme: High Small Scale Resolution 
Numerical convergence for the XLES coupled advection scheme is uncon¬ 
ditionally obtained for increasing 3D resolution (Ales cxd). This is valid 
for both coupled schemes: RK3-RK3-CDM and CN-RK3-CDM. 

The convergence for the CN-RK3-CDM scheme with increasing small 
scale resolution (Arss —)■ cx)) and two choices of constant 3D resolution 
A^les = {16,64} is additionally investigated. This kind of numerical conver¬ 
gence is not required to ensure a well dehned and converging XLES approach. 
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Figure C.6: A linear wave is propagated with constant velocity Uj = Cj over 
10 wavelengths. The coupled numerical schemes CN-RK3-CDM (dashed) 
and RK3-RK3-CDM (solid) with N = A^rss = A^les = {32, 64} are shown. 
The black line is the analytical result. 


A phase shift <F describes the numerical dispersion and thus measures the 
numerical accuracy (see hgure C.7). 

We hnd the numerical accuracy of the coupled advection scheme increas¬ 
ing up to a resolution ratio A^rss/A^les ^ 10. With higher ratios the numer¬ 
ical error of the coarse resolved scheme (RK3-CDM) dominates the coupled 
numerical error. 

Thus the numerical accuracy is increased by increasing the resolved small 
scales within XLES up to some limit in addition to possibly resolving addi¬ 
tional physical effects. 


Appendix C.3. Coupled Advection Scheme: Spectral Scale Separation 

In this section the CN-RK3-CDM scheme demonstrates the preservation 
of (spectrally separated) large and small scale properties which are repre¬ 
sented w.l.o.g. in XLES-grid 2. 

A multi-scale wave with the initial condition 


■U 2 ,i = sin(47ra;2,i) + 0.2 sin(1287ra;2,i) 
“3,1 = [K^][h]u2p = sin(47ra;3R) 
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(C.3) 

(C.4) 





Figure C.7: $ is the phase difference between the analytical solution and the 
numerical result. Case 1: RK3-RK3-CDM scheme with N = Ni^es = -^rss- 
Case 2: CN-RK3-CDM scheme with N = N^ss and A^les = {16,64}. For 


$ ~ iV 


-2 


a CN-CDM scheme Tsai et al. (2002) analytically predict the behavior 
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Figure C.8: A multi-scale wave with the initial condition in Eq. (C.3) (lines) 
propagates for 10 wavelengths. The full domain (small box) and one wave 
peak are shown. The big points corresponds to the 3D large scale held 
(A^les = 64 cells). The small points corresponds to the resolved small scales 
(A^rss = 4096). The black line corresponds to the analytical result. 


is propagated for 10 wavelengths using A^les = 64 and A^rss = 4096 cells in 
the X 2 -direction (advanced by Eq. (45) and (44)). 

Note that the higher mode 0.2 sin(1287rx2,i) is only resolved within XLES- 
grid 2 while XLES-grid 3 only captures the lower mode. 

The box hltered and small scale properties are preserved, but little dis¬ 
persed due to numerical effects (see hgure C.8). 
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